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Abstract 

We present the pr spectra, elliptic flow (^2) and coalescence parameters B-2, for d, 
d (1 < pT < 4 GeV/c) and B3 for ^He, ^He (2 < pr < 6 GeV/c) produced at 
mid-rapidity in Au+Au collisions at ^/s^ = 200 GeV. The results arc measured 
in the STAR experiment at RHIC. The spectra of the light nuclei show softer px 
distributions than calculations from a blast-wave model in which the parameters 
were fixed from pion, kaon and proton V2{pt) and pr distributions. The coalescence 
volume is found to track with pion HBT results for different collision geometries. 
The V2 measurement for d{d) as a function of transverse momentum pT is found 
to follow an approximate atomic mass number (A) scaling while that of ^Hei^He) 
deviates slightly from the scaling. A negative V2 has been observed for d at low pT, 
consistent with large radial flow in Au+Au collisions. 



Key words: 
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1 Introduction 



In relativistic heavy ion collisions, light nuclei and anti-nuclei are formed 
through coalescence of produced nucleons and anti-nucleons or participant 
nucleons [T|2|3] . Since the binding energy is small, this formation process can 
only happen at a late stage of the evolution of the system when interactions be- 
tween nucleons and other particles are weak. This process is called final-state 
coalescence [T]B] . The coalescence probability is related to the local nucleon 
density. Therefore, the production of light nuclei provides a tool to measure 
collective motion and freeze-out properties, such as particle density [5] and 
correlation volume. 

Invariant yields for the production of nuclei can be related |T] to the primordial 
yields of nucleons by Equation [H 

rf^ - ^^^^'^^ ^^"^^ ^ ^^^^^^^ 



where Na, Np, and denote the yields of the particular nucleus, and of its 
constituent protons and neutrons, respectively. Ba is the coalescence parame- 
ter, is the invariant yield of nucleons or nuclei; A and Z are the atomic 
mass number and atomic number, respectively; p^, Pp and Pn are the momenta 
of the nuclei, protons and neutrons, respectively, where pa = A- pp is assumed. 
The coalescence parameter, Ba, is related to the freeze-out correlation vol- 
ume pj: Ba oc Vj~^. For an expanding fireball, the effective homogeneous 
coalescence volume decreases with transverse mass and temperature (Mf/T) 
as detailed in Ref. [2]. 

On the other hand, a blast-wave model is often used to describe the spectra 
of identified particles produced in relativistic heavy ion collisions piT] . In this 
model, the particle spectra are determined by global parameters (temperature, 
flow proflle) and particle mass. Up to now, hydrodynamic model simulations 
or blast-wave flts to the data only included elementary particles without any 
composite elements, such as nuclei |6]. Experimental measurements can be 
used to provide insights into whether the homogeneous volume (as in a coa- 
lescence picture) or mass (as in a blast-wave model) has a bigger influence on 
the yields of nuclei in heavy ion collisions. 

In the case of deuteron production at low px , B2 oc l/Vf. The ratio of deuteron 
yield over proton yield at the same pr/A (denoted as d/p, and d/p for anti- 
deuteron), is proportional to the baryon density [5|8] . This is analogous to the 
deuteron to hydrogen ratio (D/H) measurements of Big-Bang nucleosynthesis 
(BBN), which is a very sensitive probe of baryon abundance in the early uni- 
verse piTU] . Although both processes are sensitive to the baryon density, the 
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processes themselves are very different. In coalescence, a proton and a neutron 
form a deuteron due to the overlap of the quantum wave function in a dilute 
QCD medium, while in BBN the deuteron production is through p(n,7)D 
photo-production. Nevertheless, it is interesting to investigate whether the 
coalescence at much higher temperature and density, which existed also at 
the pre-BBN stage of the early universe, produces a dramatically different 
deuteron abundance when compared to the photo-production of deuterium in 
the standard BBN model. 

Coalescence has been generalized as a mechanism for partons to form hadrons 
at an early stage of heavy ion collisions PTIIT^ . It is experimentally difficult 
to study how local correlations and energy/entropy play a role in coalescence 
at the partonic level since the constituents are not directly observable. By 
studying the spectra and elliptic flow of nuclei and comparing to those of 
their constituents (nucleons), we have a better handle on how coalescence 
happens. It has been found that elliptic flow (^2) of hadrons is very sensitive 
to their underlying partonic degrees of freedom p^tfTTj . An analog in nucleus 
production is to test if nucleus elliptic flow scales with the atomic mass number 
(A). Deviations may point to a non-uniform nucleon distribution with respect 
to the collision geometry or large nuclear density fluctuations. This will provide 
valuable data for studying the freeze-out dynamics and coalescence mechanism 
in general. 

Previously, measurements of nucleus production at RHIC suffered from low 
statistics. The first observation of anti-helium [H] at RHIC and measure- 
ments of coalescence properties p!3] with (anti-)deuterons at RHIC have been 
previously published. More recently, increasing RHIC luminosity and detector 
upgrades allow more precise measurements of yield and elliptic flow of nuclei. 
In this paper, we present spectra and f 2 of the deuteron, ^He, and their anti- 
particles {d and ^He). We compare volumes extracted from nucleus production 
and from pion HBT. First measurements of ^He f 2 and high statistics d V2 at 
low pt will be presented and their implications will be discussed. 



2 Experiment and analysis 

The data presented here are obtained from the Time Projection Chamber 
(TPC) p] and the Time-Of-Fhght (TOP) detectors [17] in the STAR exper- 
iment at RHIC in the year 2004. The TOP detector system was a prototype 
module covering vr/SO rad in azimuth and — 1.0<r7<0in pseudorapid- 
ity. These analyses used a data sample of 25 million central triggered events 
(0-12% centrality) and 24 million minimum-bias triggered events (0-80% cen- 
trality) and TOP information is available for 16 million central and 15 mil- 
lion minimum-bias events. The selection of centrality is based on the Glauber 
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Fig. 

and ^He (dashed line), (c) nad distribution of d at 0.7 < pT < 1.0 GeV/c with a 
Gaussian fit plus an exponential background, (d) The distribution of squared mass 
(m^ = (p//3/7)2) for d from the TOF after TPC dE/dx selections at 2.5 < pr < 3.0 
GeV/c, with a Gaussian fit plus a linear background. 



model and is described in Ref. [H]. Measurements of the ionization energy 
loss (dE/dx) of charged tracks in the TPC gas are used to identify protons, 
deuterons, ^He and their anti-particles. By combining the particle identifica- 
tion capability of dE/dx from the TPC and velocity from the TOF, pions and 
protons can be identified in 0.3 < < 12 GeV/c pWTj . In our analysis, 
deuterons and anti-deuterons are identified by TPC for pt < ^ GeV/c and 
by TOF in the range 1 < < 4 GeV/c. However, at low px (< 1 GeV/c), 
primary deuterons are overwhelmed by background from knock-out deuterons 
from the beam pipe and inner detector material, which are difficult to sepa- 
rate from collision products. As a result, only anti-deuterons are counted as 
collision products. Because the absorption effect and detector efficiency for 
(i at pt < 1 GeV/c are not well simulated in GEANT, and other produced 
particles in the same event produce large deuteron background from detector 
material and beam pipe, only identified d in this pr range are used for ellip- 
tic flow calculation. The TPC is used to identify ^He and ^He in the range 
2<pt<Q GeV/c. 

Figure [H illustrates the particle identification techniques and methods. Panel 
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(a) shows the dE/dx of charged tracks as a function of |rigidity| {rigidity = 
momentum/ charge) measured by the TPC at —1 < 77 < 1. Panel (b) shows 
Z {Z = \og{{dE/dx)\measure/{dE/dx)\predict)) distribution for ^He and ^He 
signals, where {dE / dx)\predict is a function accounting for the curvature of 
dE/dx versus momentum p^][20] . With tight track quality selections imposed 
by requiring the number of TPC points in a track to be greater than 25 
and the distance of closest approach to the event vertex to be less than 1 
cm, the ^Hei^He) signals are essentially background free. Panel (c) shows 
mjd (standard Gaussian deviation from {dE/dx)\pr edict) distribution for d at 
0.7 < pt < 1.0 GeV/c. The signal was fit with a Gaussian function and an 
exponential background. As can be seen in panels (b) and (c) of Figure 1, the 
peak centroids in Z and nad are shifted from zero. The large shift in ^He is 
related to the higher ionization density for particles with charge number >1. 
See Ref. [2T] and Fig. 14 of Ref. [20J. It may be changed by the effects during 
ion collection. Since this offset does not affect the extracted yields, no attempt 
was made to more carefully match the overall magnitude of dE/dx. Panel (d) 
shows distribution for d at 2.5 < px < 3.0 GeV/c measured by the TOF 
after the dE/dx cut (|n(Td| < 2, in which nad denotes the standard Gaussian 
dE/dx deviation from the expected deuteron value.) [TfTO] . The signal was 
fit with a Gaussian function plus a linear background. The acceptance and 
tracking efficiencies were studied by Monte Carlo GEANT simulations of the 
STAR detector [16]. 



3 Results 

3. 1 Spectra and the coalescence parameters 

The blast-wave (BW) model has been successfully applied to reproduce the 
vr, p, and A spectra [7], and elliptic flow as well as tt HBT correlations as 
a function of transverse momentum and centrality [6j. This model calculates 
particle production properties by assuming a parameterized fireball expansion 
after the collision. The model also assumes local thermal equilibrium with an 
expansion velocity profile as a function of transverse radius, modulated by 
an azimuthal density distribution [6]. Predictions from this model for nucleus 
yields and elliptic fiow were compared to our data. Figure [2] shows the pr 
spectra of d{d) and ^He{^He). The parameters obtained from using the blast- 
wave model [6] to describe the spectra and f 2 of vr, p are used to calculate 
d{d) and ^Hei^He) spectra, with the results shown in Figure O The spectra 
and V2 of vr, p are described simultaneously by a single set of parameters, 
which is shown in Table [1] (only the 0-80% centrality parameters are shown in 
the table). In the BW model, the nuclei are simply treated as heavier particles 
emitted from the fireball. The parameters used in the BW description are 
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Fig. 2. The left (right) plot shows d and d {^He and ^He) spectra, with a comparison 
to the blast-wave model calculations. In each plot, the upper panel shows the pT 
spectra, with the solid symbols and open symbols representing the particles and 
anti-particles, respectively. The corresponding blast-wave calculations are shown by 
solid and dashed lines. The lower panel in each plot shows the data divided by 
the blast-wave calculation. The bands show the same ratio for protons. Errors are 
statistical only. 

T, Po, P2, Rx/Ry, To and At, where T is the temperature. The freeze-out 
distribution is assumed to be infinite along the beam direction [z direction) 
and elliptical in the transverse direction {x-y plane), Rx and Ry giving the radii 
of the ellipse. The parameters po and p2 are the zeroth and second order factor 
of the flow boost along the direction perpendicular to the transverse ellipse, 
respectively. The source is assumed to emit particles over a finite duration 
in longitudinal freeze-out proper time (r = y/t'^ — z"^) peaked at Tq with a 
Gaussian distribution width Ar. A more detailed description and definition 
of the parameters and the BW formulae are given in Ref ^6j. 

In general, the spectra are similar to the model calculations. The ratio be- 
tween data and model is shown in the bottom panels of Figure [2l The model 
calculation tends to give harder spectra than seen in the data. This indicates 
that the coalescence process results in nuclei pr spectra which do not follow 
the simple mass dependence expected in the BW model, implying a deviation 
from simple thermal production. 

Table 1 



BW parameters: see text for details. 



r(MeV) 


Po 


P2 


Rx/Ry 


ro(fm/c) 


AT(fm/c) 


124.2 


0.88 


0.061 


0.89 


9.2 


0.03 



10 



n 

'o 

CM 
> 

010"^ 

I " 
OQ 



CM 

OQ 



10' 



□ 
n 

u 



► d 



-e- 



nuclei: ~ 0.8 GeV/c per nucleon 



□ 



n 



□ 



^ o 



□ 
o 



pion HBT p =0.5GeV/c 
pion HBT p^=0.2GeV/c 



° ^ □ 



O 



O 



10-' E 



10' 



50 100 150 200 250 300 350 400 



<'^art> 



Fig. 3. i?2 and yiTj together with vr HBT volume as a function of colhsion cen- 
trahty {{Npart)) in Au+Au colhsions. HBT volume is calculated from the HBT cor- 
relation lengths along the longitudinal and transverse directions. The gray bands 
represent systematic errors and the brackets show the uncertainties from the feed- 
down estimation. 



Figure [3] shows the coalescence parameters B2 and y/Rs, which are derived 
from p(j)), d{d) and ^He{^He) spectra by Equation [H Here the proton and 
anti-proton spectra are taken from Ref. [7]. The spectra have been cor- 
rected for feed-down from A(A) and weak decays |i7j. As mentioned before, 
Ba oc 1//-^, so B2 oc l/Vf and ^3 oc l/Vf. Therefore B2 for d{d) should be 
proportional to \/B^i for ^Hei^He) if the correlation volumes for d{d) and 
^Hei^He) are similar. Both B2 and -B3 show strong centrality dependence. In 
more central collisions, the smaller coalescence parameter indicates that the 
correlation volume at thermal freeze-out is larger than for peripheral collisions. 
This also means that the correlation length in nucleus coalescence grows with 
the system size. 

Figure [3] also shows a comparison with the results of pion HBT data. To 
calculate the freeze-out volume from HBT measurements, we use the following 
expression: 

Vf = (27r)3/2 X X i?Le 

where Vf is the freeze-out volume [2] and Riong and Rside are the longitudi- 
nal and sideward radii, respectively, assuming a density distribution of Gaus- 
sian shape in all three dimensions. The Riong and Rside values are taken from 
Ref. [6|22] [kx = 0.2 GeV/c). The d{d) and ^He{^He) transverse momentum 
ranges are 1.5 < px < 2.0 GeV/c and 2.0 < px < 2.5 GeV/c, respectively. The 
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HBT data is chosen at the closest px to the pr/A for the nuclei coalescence 
data throughout all of the centrality bins. The observations that the B2 and 



^/B^ coalescence parameters are proportional to 1/Vf from pion HBT over the 

full range of centrality considered indicates that the freeze-out volume for the 
nuclei is closely related to that for pions. 

In the coalescence model |2], the proportionality factors quantitatively con- 
necting the B2 and parameters to the homogeneous volume depend on 
flow profile, temperature, correction factors due to quantum wave functions, 
and other detailed assumptions of the coalescence models. A precise extrac- 
tion of these mo del- dependent factors from data will be possible in the future 
when the large uncertainty on B2 and S3 is reduced with the improvement of 
weak-decay feed-down correction to the proton yields [7|. 



3.2 Elliptic flow parameter V2 

The elliptic flow parameter, V2, is the second order Fourier coefficient of the 
azimuthal distribution of the produced nuclei relative to the reaction plane 
of the initial nucleus-nucleus collision. The event-plane method was used to 
obtain the V2 of nuclei [23j, with the event plane resolution used for correc- 
tion (calculated using the sub-event method p3]) of 76% for minimum bias 
triggered events and 68% for central triggered events. In the following discus- 
sions only statistical errors will be shown since the low statistics for deuteron 
and helium nuclei result in these errors being much larger than the estimated 
systematic errors. 

Figure H] panel (a) shows V2 as a function of pt for d + d, ^He + ^He and 
d in minimum-bias collisions. The results with both V2 and pt scaled by A 
are shown in panel (b). As mentioned in the previous subsection, the vr, K, 
p spectra and V2 are all described by a single set of blast-wave parameters 
and then the spectra and V2 of d{d), ^He{^He) are calculated. The blast- wave 
results for the deuteron i^He) V2 are shown as the solid (dashed) line. As a 
comparison, the p and A + A f 2 [22] are superimposed on the plot. Results for 
deuteron V2 from coalescence and dynamic models [12] are also shown. The 
data suggest that the d + d and baryon f 2 seem to follow the A scaling within 
errors, indicating that the d+d are formed through the coalescence of p(p) and 
n(n) just before thermal freeze-out. However, the scaled ^He + ^He V2 appears 
to deviate a bit more from the blast- wave calculated f 2 . To quantify the degree 
of agreement with the scaling, we performed a analysis by comparing the 
nucleus V2 to the curves which describe the baryon data. The per degree of 
freedom for deuteron V2 is 3.1/2 while that of ^He is 4.1/2. There appears to 
be a slightly larger deviation for heavier elements. 
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Fig. 4. (a) The elliptic flow parameter V2 from minimum bias collisions as a function 
of pt for ^He + ^He (triangles), d + d (filled circles), and d (open circles); the solid 
(dashed) line represents the deuteron {^He) V2 calculated by the blast- wave model, 
(b) d + d and ^He + ^He V2 as a function of pT, both V2 and pT have been scaled by 
A. Errors are statistical only, p (open diamonds) and A+A (solid stars) V2 are shown 
in the plot as a comparison. Coalescence and dynamic simulation for deuteron from 
Ref. [12] are also shown, (c) Low pT d V2/A (open circles) as a function of centrality 
fraction (0-10%, 10-20%, 20-40%, 40-80%, respectively). Errors are statistical 
only, p V2 is also shown as open diamonds. Blast-wave calculations are show as solid 
(d) and dashed lines (p).The 2 subpanels are for different px ranges. 

The anti-deuteron V2/A a.s a function of centrality fraction is shown in Figure H] 
panel (c). The upper and lower sub-panels represent results for two different 
regions of pt- The d is observed to have a negative f 2 in central and mid-central 
collisions in the transverse momentum range of 0.2 < px < 0.7 GeV/c. This 
negative V2 is consistent with a large radial flow, as the blast-wave calculations 
show. At the same pr/A where the d is negative, the p V2 is consistent both 
with zero and with the d V2, due to the large uncertainties. The blast- wave 
parameters published in Ref. [22] are used to calculate p and d flow and the 
calculated results are shown in the figure. The pt ranges of the p data points 
are selected to match approximately the same pr/A center of the d points used 
here. Though the blast-wave model predicts the generic feature of negative 
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V2, quantitative agreement between data and model throughout the entire 
centrahty and pt range is lacking. 

3.3 Baryon density at fiB = 

As discussed in the introductory section, the d/p ratio is proportional to 
baryon density. In the collider configuration, the production of low-p^- deuterons 
is often overwhelmed by background deuterons from the interaction of ener- 
getic hadrons (pions and protons) with detector material close to the beam 
(beam pipe, etc.). When the net baryon density is close to zero, anti-deuterons 
can be used as a measure of deuteron production. Reference [5J shows that the 
baryon density in 7P, pp, pA and AA collisions follows a universal distribution 
as a function of beam energy and can be described by statistical processes. 
At zero baryon chemical potential (/i^ = 0), the d/p ratio and the d/p ratio 
are identical and the measurements from all systems considered are consistent 
with each other. For the 5 data points, which are closest to the /i^ = con- 
dition (Au+Au collisions at ^/s^ = 200 GeV (1 point from STAR and the 
other from PHENIX), e+e" ^ ggg at = 10 GeV, jp at y/s = 200 GeV, 
and p + p at y/s = 1.8 TeV), the average value is d/p = (7.6 ± 0.8) x 10""^. 
We note that the D/H value of 2.8 ± 0.2) x 10~^ obtained from Big Bang 
nucleosynthesis in the evolution of the Universe [9l|T0] is about 4% of what is 
obtained in higher energy processes. 



4 Summary 

Using the particle identification capabilities of the STAR TPC and TOF de- 
tectors, we have measured the d{d) and ^He{^He) pr spectra. The extracted 
coalescence parameters B2 and \/B^ have similar values. They have smaller 
values for more central collisions, which is consistent with an increasing source 
size with an increase in collision centrality. The nuclei coalescence parameters 
are proportional to the inverse of the freeze-out volume estimated using the 
pion HBT radii in all centrality classes. The spectra of nuclei with A = 2 and 
3 are in general described by the blast-wave model, which was used to describe 
the spectrum, elliptic fiow, and HBT results. However, the model overpredicts 
the radial fiow, which implies that the coalescence process is different from 
the simple mass effect as assumed in the blast-wave model. 

We have also measured the V2 of d{d) and ^He{^He). The V2 values of d{d) 
when scaled by atomic mass number A, follows the baryon V2, thereby provid- 
ing evidence of d{d) formation through final-state coalescence of nucleons. We 
also observed the V2 for d to be negative at low pr in the mid-central collisions. 
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Comparison with blast-wave calculations shows this is consistent with a large 
radial flow in Au+Au collisions at a small impact parameter. 
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